Quantum process tomography of a single sòlid state qubit 
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We present an example of quantum process tomography performed on a single sòlid state qubit. 
The qubit used is two energy levels of the triplet state in the Nitrogen- Vacancy defect in Diamond. 
Quantum process tomography is applied to a qubit which has been allowed to decohere for three 
different time periods. In each case the process is found in terms of the % matrix representat ion and 
the afïine map representat ion. The discrepancy between experiment ally estimated process and the 
closest physically vàlid process is noted. 

PACS numbers: 03.67.-a, 61.48.+C, 73.21.-b 



Quantum Process Tomography (QPT) [J, Q is a 
method designed to experimentally determine complete 
information about a quantum process (quantum chan- 
nel). This knowledge can aid in identifying any sources 
of errors encountered in performing a quantum gate Q. 
QPT has been performed in nuclear magnètic resonance 
in optical systems 0, Q, and, to the best of our 
knowledge, we here present the first such demonstration 
of QPT for an individual solid-state based qubit. Meth- 
ods have been developed to determine a process from a 
tomographically incomplete set of measurements as 
well as techniques to ascertain the master equation de- 
scribing time evolution of the system 0, . There are 
two different methods with which one can perform QPT. 
The first method is to examine the effect of the unknown 
operation on a complete basis of input states. The sec- 
ond method is to examine the effect of the channel on 
a larger state composed of the original system and an 
ancilla system 0, Q • In either case state tomography is 
used to ascertain the output state (s) and this informa- 
tion is used in conjunction with the known input state (s) 
to derive the process. The former method will be applied 
in what follows. 

The Nitrogen- Vacancy defect consists of a substitu- 
tional nitrogen impurity next to a vacancy in the dia- 
mond lattice. The resulting electron configuration pro- 
duces a spin triplet (S=l) state which can be manipu- 
lated by electron spin resonance and read out by opti- 
cally detected magnètic resonance 

(Fig. l(a)). At low temperature, the spin longitudinal 
relaxat ion time, Ti, is on the order of seconds and thus a 
single electron spin state can be detected . Transverse 
relaxation times (or decoherence times), T2, of up to 60 
microseconds have been reported in literature |l4|, for 
samples with low nitrogen concentration. A controlled 
two-qubit quantum gate in which the vacancy electron 
spin is hyperfine coupled to a nearby C 13 nucleus, has 
recently been performed using this system 0. 

Using quantum process tomography one can determine 
a completely positive and trace-preserving map £ which 



represents the process acting on an arbitrary input state 
P : 
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where the À m are a chosen basis for operators acting on 
p, and d = 2 n for an n qubit system. The matrix \ 
completely describes the process £ , and can be recon- 
structed from experimental quantum state tomographic 
measurements. To obtain x it is first necessary to ap- 
ply the process to each member of a complete basis of 
input states (e.g p s = |^><^U^> = {|0>, |1>, ^-(|0) + 
|l)),^(|0)+i|l))}, for a single qubit) and perform state 
tomography on the resulting output states. From these 
measurements Xjk can be specified m the relation 
£{pj) = XjkPk where pj form a basis for density matri- 
ces. In order to determine \ from À one operates on A 
with the pseudoinverse of /?, where /3 is derived theoret- 
ically from the relation AmPjA^ — fljj^pk- Using this 
last relation and (QJ we can see À^- = £ mn (3\ 
so inverting (3 gives us as required. 

The QPT experiment was performed at room temper- 
ature using diamond nanocrytsals obtained from type Ib 
synthetic diamond. Coherence times on the order of a mi- 
crosecond have been reported previously for this type of 
diamond. In order to perform many repetitions in a rea- 
sonable amount of time, a sample with a relat ively short 
coherence time was chosen. Diamond nanocrystals were 
spin coated on a glass substrate, and single nanocrystals 
were observed with a homebuilt sample-scanning confo- 
cal microscope. In order to ensure the presence of a sin- 
gle defect in the làser focus, the second-order coherence 
was measured using Hanbury-Brown and Twiss interfer- 
ometer and the contrast of the antibunching depth was 
determined. Microwaves were coupled to the sample by a 
ESR microresonator connected to a 40W travelling wave 
tube amplifier. 

The bàsic energy level scheme of the molecule relevant 
to the experiment is shown in Fig. l(a). The ground 
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electrònic triplet state consists of the three spin sublevels 
corresponding to spin projections ms = ±1,0. Sublevels 
with ms = and ms = 1, were selected as qubit states. 



Làser 



form: £4 = 



1 

t E 



where E is a real 3x3 matrix 
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FIG. 1: (a) Scheme of energy levels of NV center in weak 
magnètic field. Because of the C3v symmetry of the NV de- 
fect, the ms — ±1 sublevels remain degenerate in zero field 
and are split in non-zero magnètic field. (b) Experimental 
pulse sequence. Note that the làser was switched off during 
spin manipulation in order to avoid decoherence related to 
continuous spin measurement via optical excitation. 



In a preliminary step of the experiment, the times re- 
quired to perform ir and Ç pulses were determined using 
Rabi oscillations. After that the experimental sequence 
shown in Fig. l(b) was applied. As a first step the qubit 
was initialized in the ms = state. This was achieved 
by optical pumping, which results in a strong spin po- 
larization, corresponding to at least 70% population of 
the ms = sublevel. Optically induced polarization is 
related to spin-selective intersystem crossing from pho- 
toexcited 3 E triplet state to metastable singlet state 1 A 
A complete basis of four input states was then prepared 
using microwave pulses resonant with qubit trànsit ions. 
The po state was obtained directly by optical pumping 
and three remaining input states pi, p2 and ps were ob- 
tained by application of suitable tt or Ç pulses. Each 
of these states was left to decohere for three different 
time intervals: 20 ns, 40 ns and 80 ns. As a last step, 
measurements of the diagonal and off-diagonal elements 
of the density matrix were performed. Estimates of the 
density matrix elements were extracted from experimen- 
tal data (Rabi oscillations) using the maximum entropy 
(MaxEnt) technique This method returns a phys- 

ically vàlid density matrix which satisfies, as closely as 
possible, the expectation vàlues of measured observables. 
In cases where only incomplete knowledge of the output 
state is known, an additional constraint is used; the re- 
constructed state must also have the maximum allowable 
von Neumann entropy. 

The state of an arbitrary qubit can be expressed as 
p \{X + r • a) <-> \(l,x,y,z) where 3 {cr x ,a y ,a z } 
and x 2 ± y 2 ± z 2 < 1. In this basis an arbitrary trace- 
preserving process £ can be written in the afflne map 



(responsible for deformation and rotation of the Bloch 
sphere) and t denotes displacement from f = (0,0,0). 

Using this picture of quantum processes , one can 
reconstruct the map £ by examining the transformations 
undergone by a complete basis of states: po = \l, p x = 
\{1 ± a x ), p y = \{1 ± (j y ) and p z = \{l+(j z ). For 

brevity we write £U(po) = m 'i £ A(Px) = x',£ A {p y ) y' 
and £a(Pz) = z' • The transformation £ can then be 
expressed in the afflne form by 
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m x' — m y 



A map is completely positive 0, and trace preserv- 
ing if and only if it can be written in the Kraus (op- 
erator sum) representation: £(p) = Àj^pÀ^. where 
^ k À\Àk = 1>. Physically, the complete positivity re- 
quirement states that, for a one-qubit process £ and an 
entangled Bell state |^ + ) = ^(|00) ± |11)), {£ 
must also be a vàlid two qubit state. 

The process obtained directly from experimental data 
is often unphysical i.e. non-trace-preserving or not com- 
pletely positive. In such cases it is necessary to search for 
a physical process which is closest in some sense to the 
measured experimental results. In this case we used a 
least squares fit between the experimentally determined 
X and a Hermitian parametrization ^| of a physical \ 
while enforcing complete positivity and trace preserva- 
tion. 

The graphical representation of a process as a defor- 
mation and rotation of the bloch sphere (Fig. 2-4) can 
give crude but obvious insights into the physicality or 
otherwise of the estimated experimental process (protru- 
sions outside the unit bloch sphere as well as ellipsoids 
which are "flattened" in one dimension betray a violation 
of the trace preservation and complete positivity condi- 
tions) and also into the fidelity of the desired unitary pro- 
cess. It is important to note that not all ellipsoids within 
the bloch sphere can be obtained by a completely posi- 
tive map. Transverse relaxation (the dephasing channel) 
is depicted as the Bloch sphere collapsing to the z axis. 
As expected this is seen to be dominant decoherence pro- 
cess. Longitudinal relaxation (amplitude damping) is not 
discernible on the timescales used in the experiment. 

It is possible to convert the process £ into a de nsity 
matrix, pg, via the Jamiolkowski isomorphism 0, liof : 
P£ = [J®£](|*)(*|) where |$) = £\ \j)\j)/Vd and 
{\j}} is an orthonormal basis. When the process £ is 
physical, one then obtains a physically vàlid pg which 
can be compared to the ideal process (converted to pid) 
using distance measures on quantum states 0. The 
trace distance between density matrices pid and pg is 
D(p id ,ps) = \tr | pid - PE |, where | X |= VxVC. 
Similarly one can define the Fidelity: F(p i d 1 pg) = 
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Decoherence 
time 
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Dp ro 


20ns 


0.101 


0.050 


0.066 


0.056 


40ns 


0.110 


0.059 


0.075 


0.062 


80ns 


0.175 


0.075 


0.110 


0.096 



TABLE I: The disparity between the estimated process from 
the experimental measurements (x) and the "nearest" physi- 
cally vàlid process (%). 
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ír (y^y/pïd Pe^fpïd) and, using this definition, 
Bures Mètric B(p id ,pg) = ^2- 2^F(p idl p s ) and 

C mètric C(p id ,pg) = a/1 - F(p id ,pe). 

An unphysical process, however, can lead to an unphys- 
ical ps, possibly resulting in a process fidelity which 
is greater than one. The application of the preceding 
fidelity-based distance measures can, therefore, produce 
nonsensical results. In such cases it is necessary to use 
other techniques in order to estimate the disparity be- 
tween x an d x f° r example. If one defines X = x~ X? then 
possible measures are the matrix p-norms (p = l,2,oo) 
of X and the Frobenius norm of X (||X||i? ro ) as well as 
the trace distance (D pro ). 

As stated above, these quant ities gives a measure of 
how well x describes the experimental results. The fi- 
delities corresponding to these measures for the different 
decoherence durations are presented in Table I (note that 
||X||i = || X Hoc because X is Hermitian and that the ba- 
sis operators used for x (i- e - the operation elements) are 

{(J0,(Tx, —ÍO'yiO'z})' 

We have presented the first quantum process tomo- 
graphic analysis of an individual single sòlid- state qubit, 
a Nitrogen Vacancy centre in Diamond. This analysis 
is only possible due to the enormous advances made in 
recent years in single- molecule spectroscopy 0, HM , 
where the resultant ODMR technique provides us here 
with high- fidelity single-qubit readout. As quantum de- 
vices develop and increase in size, the task of "debug- 
ging" the device, or actively identifying the noise present 
in the device, will pose significant challenges. The work 
presented here represents an initial step towards the test- 
ing of quantum devices in solid-state. 

M.H. acknowledges funding from the Irish Research 
Council for Science, Engineering and Technology "Em- 
bark" Initiative. The work was also supported by the 
EC FP5 QIPC project QIPDDF-ROSES. 



[1] M.A. Nielsen and I.L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 



FIG. 2: 20ns results: (a) Bloch sphere representation of process 
from experimental data. (b) Bloch sphere representation of 
process obtained from x- 




FIG. 3: 40ns results: (a) Bloch sphere representation of process 
from experimental data. (b) Bloch sphere representation of 
process obtained from x- 
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process obtained from %. 
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